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We question whether the measurement based quantum computing algorithm is in fact Grover's 
algorithm or simply a similar oracular search method. The two algorithms share several qualitative 
features especially in the case of the trivial 4 element search, which is the largest size photonic search 
algorithm that has been experimentally implemented to date. This has led some to refer to both 
substantiations as Grover's algorithm. We compare multiple features of the two algorithms including 
the behavior of the oracle tags and the entanglement dynamics, both qualitatively and quantitatively. 
We find significant and fundamental differences in the operation of the two algorithms, particularly 
in cases involving searches on more than four elements. 
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I. INTRODUCTION 



L. Grover first showed that a speed up over classical database searches was possible using quantum mechanics. 
His search algorithm, now known as Grover's algorithm [1], along with Shor's factoring algorithm [2] helped to spur 
I the interest in all manners of quantum computing. Multiple implementations of universal quantum computing have 
been found and tested, such as Linear Optical Quantum Computing LOQC 3'], Adiabatic Quantum Computing 
(AQC) [4 hand Measurement Based Quantum Computing (MBQC) or equivalently cluster state quantum computing 
(CSQC) H @ just to name a few. Grover's algorithm has been shown by C. Zalka Q to be the optimal means of 
searching an unsorted list via amplitude amplification and, with some minor modifications, has been shown to have 
near perfect fidelity in the limit of large system size N Q. We will compare and contrast the MBQC search with 
Grover's algorithm in the photonic circuit model. In Section II we describe the initial assumption that led to this line 
, of inquiry. In Section III we describe in detail the MBQC implementation and our numerical framework. Section IV 
explicitly defines the two algorithms. Section V describes one of the main differences between the algorithms, namely 
| the tagging operation in the MBQC search. Section VI deals with amplitude amplification and compares the speed 
of the two algorithms. Section VII discusses the differences in the entanglement dynamics. Section VIII asks several 
unresolved questions about the MBQC search. Section IX summarizes our results and makes concluding remarks. 
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II. ASSUMPTION THAT CIRCUIT MODEL AND MEASUREMENT BASED MODEL GROVER 

SEARCH ALGORITHMS ARE THE SAME. 



In almost every implementation of quantum computing one of the first algorithms that people have developed and 
tested is an unstructured quantum search algorithm. In the photonic circuit model of quantum computing [9] this 
algorithm is structurally very similar to that first proposed by Grover and as such is called Grover's algorithm. Since 
the development of MBQC, and particularly in the photonic implementation of MBQC, also called the cluster state 
model, parallels have be drawn between the photonic circuit model and MBQC model. This includes comparing the 
circuit model of Grover's algorithm to a cluster state model of a 4 element quantum search [n| ■ Given the difficulty of 
creating large cluster states this four element search is the largest that has been preformed to date. The, 2™ = 2 2 = 4 
element search requires n 2 = 2 2 = 4 qubits arranged in a square or "box" cluster state [To| . The next larger 2 3 = 8 
element search requires a significantly more difficult to produce 3x3 square cluster state and so on. 

Using the 4 qubit box state, the MBQC search appears to be very similar to that of Grover. The input state is 
similar, the desired output state is similar, the number of iterations is the same and the oracle's tagging operation 
appears similar. This has lead many people to call the MBQC quantum search, Grover's algorithm. However there is 
often a small caveat that is overlooked. As Zeilinger et al. wrote "Remarkably, the inversion-about-the-mean process 
is 'hard- wired' into the structure of the cluster state and is automatically implemented" [l(| ■ We found this automatic 
implementation of an inversion about a mean, or equivalently amplitude amplification (when combined with tagging), 
to be a imprecise claim and as such have investigated the measurement based quantum search algorithm numerically 
for larger systems. We have found that a significant number of non-trivial differences exist between the standard 
description of Grover's algorithm and its implementation in the circuit model and the MBQC search algorithm for 



2 



larger systems. This leads us to question whether the MBQC search algorithm is some variant of Grover's algorithm, 
or rather a different search algorithm altogether. 



III. THE MEASUREMENT BASED QUANTUM COMPUTING AND THE NUMERICAL FORMALISM 

Measurement based quantum computing (MBQC) operates by first generating highly entangled multi-qubit states 
called cluster states. Computations are then preformed by a time ordered sequence of measurements on individual 
qubits in this entangled resource. Some cluster states, such as the 2 dimensional square cluster state (see Fig. ([1])), 
have been proven to be able to perform universal computation assuming the initial cluster state is large enough. 

To begin the construction of the cluster state, we first form the Kronecker product of N photons each initialized to 
the |+) = (|0) + | l))/v / 2 state. 
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However, to complete the construction of the cluster state we must create the entanglement resource between the 
qubits by effectively performing the CZ operation between nearest neighbors. This entangling operation between two 
qubits can equivalently come from other sources such as spontaneous parametric down conversion (SPDC) llj. The 
four CZ operations for the box cluster state (see Fig.(QJ,)) are shown in Eq.@, 
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>) = l^bo.) = H CZiMimt) = CZ 1 , 2 CZ lt3 CZ 2A CZ 3A \i(>), 



(2) 



where CZ%j is a 4 x 4 matrix with entries (— l) 6 ' bj for the binary computational basis entries fcj bj € {0, 1}. Here we 
note that i,j are the labels of the control and target qubits respectively in the CZ operation, nn is a list of all nearest 
neighbors and the that order of the CZ operations is irrelevant. Two of the cluster states of interest in the remainder 
of this paper are then shown in Fig. (p}; the 4 qubit "box" state Y&box) and its next logical extension, the 9 qubit 



square cluster state \^c9l 
N = n 2 total qubits. 



The CZi j operations create the n x n square cluster state the MBQC search requires with 
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FIG. 1: Square Cluster-states. Two cluster-states able to perform an unsorted search for a) the (n = 2) 2 2 = 4 element 
search on \^box) and b) the (n = 3) 2 3 = 8 element search on \$!c9}- Each circle is a single qubit initialized in the |+) state. 
The solid lines connecting the circles indicate the action of a CZ gate between these qubits, i.e. entanglement between the 
qubits. Logical qubits are arranged in rows of the grid, initially: (1,3) and (2,4) for \9box), and (1,4,7), (2,5,8), and (3,6,9) 
for l^og). a) shows the trivial MBQC measurements using one application of an iterator I\ (large red dashes acting on columns 
of qubits) on the 4 qubit box state, b) Shows the two iterations for the 8 element search 7i,/2. All of the measurements in 
an iterator Ij can be implemented simultaneously. The outcomes of each Ij determine the basis of the next iteration Ij+i. 
(Equivalently Ij determines the Pauli correction factors that can be applied latter in post-processing) . The output state (small 
green dashes) holds the logical state at the end of the algorithm. The effective logical qubits after the first iteration I\ are (3) 
and (4) (the output states) for \^ box ), and (4,7), (5,8), and (6,9) for |*C9}- 



An algorithm then proceeds as a series of single-qubit measurements. The bases for the measurements performed at 
the (k + l)th iteration are determined by the outcomes of the previous measurements performed at iteration k. This 
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can be considered as an ordered series of Kronecker products of single qubit rotations. Equivalently we may rotate 
the measurement basis. Such a Kronecker product of rotations produces a unitary matrix. For the MBQC search 
algorithm we need only the set of rotations of s = {0, tt} [lfj. We are measuring in the basis |± s ) = (|0) ± e ls \l))/y/2, 
i.e. either the standard |±) basis (s = 0) or the |=f) basis (s — tt). 

The MBQC search proceeds by implementing measurements on the (n — 1) vertical columns of qubits, each of size 
n on the cluster state. It proceeds sequentially from left to right, leaving the last (i.e. rightmost) column of n qubits 
as the output state of the computation, see Fig. (TTJ). Operationally, each measurement B(a) is implemented by the 
product of a rotation matrix R(a) times the Hadamard matrix H , 

B(a) = R(a) H, R{a) = ( CO&a 'l V H = ± ( \ \) , (3) 

v 1 \ i i \ zsma/2 cosa/2 J ' ^fo \ 1 — 1 I v ' 

followed by a projection P onto |± s ). Note that the rightmost H in B(a) rotates the computational basis states 
{|0), |1)} (eigenstates of Z) into the states {|+), |— )} (eigenstates of X) respectively, while R(a) effects the rotations 
R(a) |±) = e ±ra / 2 |±) . In this formalism the unitary matrices that implement the measurement based search algorithm 
takes the following forms for the two search sizes of interest shown in Fig.(fT]): 

U box = B 1 {a)®B 2 ((3)®H 3 (g>H 4 , (4a) 
U C 9 = Bx{a) ® B 2 (fi) <g> B 3 (j) <g> S 4 (o) <8 B 5 (b) ® B 6 (c) ® H 7 ® H s ® Hg. (4b) 

The output state of the calculation, prior to the the projection operations, is then, 

\Vout)=U\% n ). (5) 

The output state \^ ou t) is a vector of size 2 N . The output state can divided into 2 N ~ n sections of size 2 n , each of 
which correspond to a set of specific measurements outcomes on the n — 1 columns of n qubits. In the first section of 
2™ elements of the state \^ ut), each of the prior n(n— 1) qubits have been projected "correctly" into the state |+ s ) (i.e 
|+) for s = 0, and |— ) for s = it ) for their given value of s (i.e. |* out ) = |+ ai )i® ■ • • |+ Sjv _„)w- n ® |Vw)(jv-n)+i,...,iv 
where \ipout) lN-n)+i,...,N is the output state on the last column of n qubits) yielding the output for the MBQC search 
algorithm. Each of the remaining segments, each of size 2™, correspond to "incorrect results" in which one or more 
of the prior n(n — 1) qubits has been projected into the |— s ) for their given value of s. However, like any MBQC 
algorithm these errors are known, since we know the outcomes of the measurements, and they can be corrected 
during or after the calculation with polynomial resource overhead. For full details see [12| or other texts on MBQC 
[l3T |. The essential point is that these "incorrect outputs" are calculable and can be trivially mapped back into the 
correct output based on the measurement outcomes for any size system. Therefore, we can assume without loss of 
generality that the single qubits measurements all give the "correct" output. This eliminates the need to deal with 
trivial feed forward corrections at the end of the algorithm. We can apply the projection matrix P = |+ s )(+ s | to 
each measurement (See Appendix [XI for details). The measurement operations for the 4 qubit box Eq. (|4"aj) are then, 

U box = PB^a) ® PB 2 (fi) ®H 3 ®H Al P= \+ s ){+ s \ = ( J J ) • (6) 

This projective operation populates the unitary matrix U bo x with the non-zero values in the first 2 n rows and zeros in 
the rest of the matrix. This is convenient for calculations involving large cluster states for which using dense matrices 
of size 2" x 2™ is computational expensive. We can simply assume the measurements in the MBQC scheme proceed 
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"correctly" and truncate the matrix to size 2™ x 2™ . Note that U is unitary in the absence of the projection matrix 
P. In the remainder of the work the projective measurements P are implied on the first n x [n — 1) qubits. 

We sample the 2 n outputs by simply rotating the n output qubits of the last column of the grid to the |±) basis. 
This is equivalent to preforming a simple Hadamard on each output qubit and measuring. Note that different possible 
measurement patterns on the cluster state are equivalent (right to left, left to right, top to bottom and bottom to 
top); one simply has different locations for the output state. 



IV. GROVER'S ALGORITHM VS. THE MBQC SEARCH ALGORITHM 



In this section we examine the details of Grover's algorithm and compare it to our implementation of the MBQC 
search algorithm. 
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A. Grover's algorithm 

Consider an unstructured set of Nq — 2™ elements X = {1, . . . , Nq} € {0, l}® n for which we desire to find a specific 
element Xq. For x g X let f(x) : {0, 1}®" — > {0, 1} be an externally supplied function (the oracle) defined such that 
f{x ) = 1 and f(x ^ x a ) = 0. Grover defined his algorithm with the following steps 

1) initialize the system of size Nq to the uniform amplitude, unbiased state \i/ji n it) = ^ J2x \ x ) 
: ! . 1, . ..,1)/ v / ]Vg (in the computational basis of states {|0), 

2) Repeat the following operations 0(^/Ng) times: 

a) Let \x) be any computational basis state: 

if f(xo) = 1 rotate the state by it, 
if f{x) = leave the state unaltered. 

b) Apply the diffusion transform D such that: 

Dij = 2/N G if i £ j k Da = -1 + 2/N G . 

3) Sample the state by projective measurement. 

It is common to refer to step "2a" as a tagging operation in which the phase of the desired state \xg) (or set of 
states, in general) is given a tt phase flip. Step "2b" is the so called inversion about the mean, which is essentially a 
rotation that performs an averaging of the state amplitudes. Both "a" and "b" can be written as unitary operators, 
which we will denote Uf and Ui nv . The iterator of Grover's algorithm is then simply the matrix product Ug = Ui nv ll /. 
Grover's algorithm begins from a uniform input state of size Nq = 2™ and preforms a sequence of repeated tag and 
inversion operations U G until the desired number of steps k ss 0(-\/Nq) have been performed. The state is then 
measured and the desired result or results will have an overwhelming likelihood of being measured. 

As is well known, the tagging operator U / requires knowledge of the solution element sought in order to be explicitly 



constructed [14j. This is why the oracle Uf is part of the "input" to the Grover search algorithm, and it has to be 
supplied externally (via an oracle) to the user performing the search. Here, Uf is the unitary implementation of 
the classical function /. In order to be unitary (and hence reversible), Uf acts on the product n-qubit basis states 
\x) and a single auxiliary qubit \y) as Uf \x) <X> \y) — \x) ® \y © f(x)), where © denotes binary (mod 2) addition. 
For \y) = | — ) = (|0) — \l))/y/2, explicit computation reveals that Uf \x) ® |— ) = (-l)f(x) \x) ® |— ). This phase 
kickback operation is often abbreviated as Uf \x) = (— l)/^) \x), though the use of the single auxiliary qubit is 
mandatory in order to ensure the unitarity of Uf. An explicit representation of tagging operator is given by [l4j 
Uf = \xo) (x \ (g) X 2 + J2 x ^x \ x ) ( x \ ® ^2) where X 2 and I 2 are the 2x2 single qubit Pauli X bit-flip operator and 
identity matrix respectively, which act on the auxiliary qubit. This latter expression demonstrates that the oracle 
must know the sought after solution xq in order to know where to place the single qubit bit-flip operator X 2 in the 
block-diagonal matrix Uf (where all the remaining 2x2 diagonal blocks for x ^ xq are filled with 12). A variation 
of Grover's algorithm that circumvents this requirement to "know" the solution state Xq in order to construct the 
tagging operator Uf is given in 15j which searches on pairs (or more) of tensor product states \x) ® \t), where say, x 
denotes a name, and t a telephone number. Here, searching on a supplied telephone number |io) in order to find an 
associated unknown name |xo) requires the construction of the tagging operator I x ® Uj, where I x acts unbiasedly 
as the identity operator on the unknown subspace of names and Uj acts in the subspace of telephone numbers. 



B. MBQC search algorithm 

In this section we describe the standard MBQC search algorithm. 

1) Initialize a set of -/V = n 2 qubits into an n x n square cluster state in which each qubit is in the plus state 
|+) = (|0) + \l))/y/2 and 2 n is the maximum desired number of elements to be searched. 

2) An oracle selects a set of measurements of size n x (n — 1), to be made in either s = or tt basis such that 
|± s ) = (|0)±e is |l))/V2 and: 

a) n — 1 sets of measurements are applied sequentially, in which n qubits are measured simultaneously, moving 
across the cluster state in a uniform direction (i.e. left to right, or top to bottom), and such that the the (k + l)th 
measurement set depends on the outcomes (successful or unsuccessful) of the fc-th measurement set, for all k = 
1, . . . , n — 1 iterations. 

3) Sample the final state by rotating the remaining n qubits into the |±) basis, followed by projective measurement. 

The most apparent difference between the definition of Grover's search algorithm and the MBQC algorithm is in the 
form and function of the iterator. In the circuit model of Grover's algorithm the iterator is a product of two unitary 
operators which repeatedly implements the same (i.e. constant between iterations) tagging operation, followed by the 
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same inversion about the mean operation. In the MBQC search algorithm this repetitive implementation of fixed, 
constant iterations is not mandatory. The size n "iterator" , which consists of the n measurement bases chosen by the 
oracle, need not, and in general will not, be a constant between iterations. 

V. TAGGING BY THE ORACLE IN THE MBQC SEARCH ALGORITHM 
A. The MBQC search algorithm with general non-constant tagging operations 

In Grover's algorithm the tagging operation, which is preformed by the oracle, is a unitary transformation which 
simply flips the phase of a desired element. As a non-trivial, illustrative example, we consider searching for a unique 
element in a eight element list. This tagging operation forms half of the Grover iterator. The tagging is one to one, 
meaning there is a single unique tag for each and every element in the search space. In order for the inversion about 
the mean to take place the oracle must continue to apply the same tagging operation at each step, i.e. iteration, of 
the search. As stated above this is different from the MBQC search. 

The MBQC search is unique only up to overall phases such as ± signs and factors of i, and the final output states can 
be degenerate. For the (n = 3) 2™ = 8 element search on the 3x3 grid of qubits, the the final possible measurement 
outcomes on the last column of 3 qubits are 8 fold degenerate. That is, there are 8 tagging operations that all point 
to the same final state (ignoring overall phases). In the special case of the (n = 2) 2™ = 4 element search implemented 
on a 2 x 2 grid of qubits (the box cluster state), there are 4 measurements and 4 unique output states, implying that 
the tagging operations are one-to-one. We believe that the one-to-one mapping between the tagging operation and 
the final output state exhibited in the special case of n = 2 has resulted in part of the confusion about the equality of 
Grover's and the MBQC search algorithms. We have found that a one-to-one property associating a tagging operation 
to a unique output state does not scale with the size of the cluster state. 

The MBQC search can be implemented successfully for any n with a constant (i.e. fixed) iterator. The (n = 3) 
2" = 8 element search on 9 qubits arranged on a 3 x 3 grid can be implemented with an iterator composed of three 
measurements a,b,c € {0, n}. This iterator would be applied twice to the two columns of 3 qubits, i.e. (a,b,c ; a,b,c) 
on qubits (1,2,3 ; 4,5,6) respectively. However, this is not a necessary condition, and the MBQC search could instead 
be preformed with an iterator that varies between measurement iterations (see Appendix [B] for the 3 and 4 qubit 
examples) . 

Interestingly, we have found that if one ignores, or corrects for overall phase, then not all of the n x (n — 1) 
measurements are strictly necessary. Note that this does not mean that some of the qubits are unnecessary, the size 
n of the qubit grid remains the same. As an example, we again use the (n — 3) 2™ — 8 element search on 9 qubits on 
a 3x3 grid. This search requires 3 x (3 — 1) = 6 measurements chosen by the oracle with the remaining 3 output qubits 
rotated to the |±) basis via an H® 3 operation. However, we find that we can preform the full 8 element search, i.e. 
produce every possible output state (ignoring the overall phase of the output) by varying only n = 3 measurements. 
Explicitly the sequence (0, 0, S3 ; 0, S5, sq) on qubits (1,2,3 ; 4,5,6), where implies a measurement projection onto the 
standard plus state |+), i.e B(0), and Si implies a measurement that can be varied, s = or s = it, by the oracle based 
on the desired output (i.e. the oracle can change the Si's to tag the desired output state, but doesn't need to change 
any of the measurements involving B(0)). It is important to note that not every combination of n — 3 measurements 
gives the full space of 2™ = 8 outputs. As an example, the measurement pattern (0, 0, ; S4, S5, sq) only results in 4 
distinct output states (up to overall phases), i.e. a 2 qubit search. The pattern of measurements (0,0, S3 ; 0,S5,S6) 
has a form similar to Eq. (l4"aj) . 

U C 9 = H X ®H 2 ® £3(7) ®H A ® B 5 (b) <g> B 6 (c) ®H 7 ®H$® H 9 . (7) 

This leads to the question of precisely what operation the "tagging" unitary is performing on the state. Above 
we define the operation of an MBQC algorithm as a unitary matrix constructed by taking the Kronecker product 
of rotation operations on the n(n — 1) measured qubits, multiplied by the Hadamard matrices on the final n output 
qubits. This leads to the Kronecker product of a series of matrices chosen from the set of or it rotations, or explicitly 
from Eq.©, 

fi(o)=^ = -L(J_ 1 1 ), & B ^ = F = ^= 2 (I "i 1 )- (8) 

The MBQC search operation U is then the Kronecker product of H's and F's. For the (n = 3) 3 x 3 = 9 qubit cluster 
state U has the following form given in Eq. (|9a|) and Eq. (|9bl) for the n(n — 1) = 6 measurements (a, f3,-f ; a, b, c), 
with the output qubits 7, 8 & 9 rotated into the |±) basis. Two specific examples of the set of 2™( n_1 ) = 64 unitary 
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matrices are shown below that produce the output states |1) and |4), or (1, 0, 0, 0, 0, 0, 0, 0) T and (0, 0, 0, 1, 0, 0, 0, 0) T 
respectively, up to an overall phase, 

Ueiementi = #1 ® #2 <8> H 3 ® Hi ® H 5 ® H 6 ® H 7 ® Hg ® H 9 , (9a) 
Ueiementi = H 1 ® H 2 ® F 3 ® H A ® F 5 ® F 6 ® H 7 ® Hg ® Hg . (9b) 

Here we pause to note something rather undesirable. If we consider the measurement of the vertical columns to 
be a "step" in the iteration (i.e. 1,2,3 or 4,5,6 from Fig. ([I] b)), then the MBQC algorithm does not have a hxed, 
constant iterator. The iterator will in general vary between steps and it is neither sequential nor random as can be 
seen from the first column of Table The various iterators must be carefully chosen by the oracle at each step. This 
appears very different from Grover's iterator. In Grover's algorithm the iterator is chosen by the oracle at the start of 
the calculation and then remains a constant regardless of the number of iteration steps. Table (HJ shows measurement 
settings for the 6 qubit measurements that result in each of the 8 possible outputs of the MBQC search. Note the 
output state is corrected for overall phase, so the tags are not unique. In addition, all measurements presented are 
presumed to have given the "correct" output (i.e. only entries in section 1 of \^!out) are populated), thus negating the 
need for trivial feed forward corrections typical in MBQC. 



TABLE I: Measurement patterns & resulting output (ignoring phase factors) 
for the (n = 3) 2 n = 8 MBQC search on a 3 x 3 grid of qubits. 



Measurement Sequence (T\ ; T2) 


Search Element i 


Output State \^ ut,i) 


(0,0,0 ; 0,0,0) 


1 


(1,0,0,0,0,0,0,0) 


(0,0,77 ; 0,7T,0) 


2 


(0,1,0,0,0,0,0,0) 


(0,0,0 ; 0,0,tt) 


3 


(0,0,1,0,0,0,0,0) 


(0,0,tt ; 0,7r,7r) 


4 


(0,0,0,1,0,0,0,0) 


(0,0,tt ; 0,0,0) 


5 


(0,0,0,0,1,0,0,0) 


(0,0,0 ; 0,tt,0) 


6 


(0,0,0,0,0,1,0,0) 


(0,0,tt ; 0,0,tt) 


7 


(0,0,0,0,0,0,1,0) 


(0,0,0 ; 0,7r,7r) 


8 


(0,0,0,0,0,0,0,1) 



B. The MBQC search algorithm with a single constant iterator for any outcome 

We have found a way of modifying the MBQC search algorithm such that the iterator is a constant not only between 
steps, but also constant for any desired output. We begin by noting that in the n — 3 case there is one example in 
which the iterator is always a constant, namely Ueiementi- In Eq. (|5a|) (see first row of Table UJ we can see that the 
two steps of the iterator are the same for Ueiementi, i.e. Hi ® Hi ® H 3 and H4 ® H 5 ® H e . To determine the effect of 
the tagging operation we compare two different U matrices, for instance Ueiementi and Ueiementi- We performed this 
comparison by the simplest expedient operation of subtraction, 

AU = U e l emen tl ~ Ueiementi = H® 9 - Hi ® H 2 ® F 3 ® H 4 ® F 5 ® F 6 ® H 7 ® H 8 ® H 9 . (10) 

Here we find an interesting form for the matrix AU. Many of the columns are identically zero while the remaining 
columns in AU have an absolute value twice that in any given element of U{. This implies that the two differently 
tagged matrices are identical except for the sign of the elements in several vertical columns. This implies that, 

Ueiementi ~ Ueiementi ?4 = H® 9 - (Hi ® H 2 ® F 3 ® H 4 ® F 5 ® F 6 ® H 7 ® H S ® Hg) T 8 = 0, (11a) 

T 4 = h®I 2 ® (~Z 3 ) ®h® (-Z 5 ) ® (-Z 6 ) ®I 7 ®I$® I 9 , (lib) 

where in Eq. (|llbj) we have used the relations F F* = I and F* H = i(—Z). The i in this relation leads to some of 
the tags being imaginary such as T4. Note that this i is simply an overall phase and without loss of generality can be 
factored out at any convenient point. For clarity we drop the i from the definition of Ti giving all real tags. In the 
above, T4 is real a diagonal (tagging) matrix with elements ±1. The l's correspond to the columns that are already 
identical and thus subtract to 0, while the -1 correspond to the remaining columns. U e iementi is identical to the 
product UeiementiTi. Since Ti is an involution TiTi — /, we can also write UeiementiT± = Ueiementi- However, the 
matrices Ti do not commute with Uj ([Ti, Uj] 7^ 0). A matrix Ti can be found for any matrix Ui, or equivalently for 
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any measurement pattern, i.e. (si, s 2 , S3 ; S4, S5, sg). The Ti matrix corresponding to the U e i emen ti transformation is 
equal to the identity matrix T± = I. 

As an illustrative example, we provide the explicit form of the T matrices in the case of the (n = 2) 2™ = 4 
element search on the box cluster state in Fig. [T^, (note: for the (n = 3) 2™ = 8 element search, the T matrices are 

of size 2" x 2™ = 512 x 512). The diagonal tagging matrices Ti for i 6 {1,2,3,4} corresponds to measurement 
iterators I, = (31,82) acting on qubits 1 and 2 given by: Zi = (0,0) — > T\ = I\ <X> h ® h ® I4, T 2 = (tt, 0) — > T2 
= (-Zi) ®I 2 <E)I 3 ® h, T 3 = (0, tt) -> T 3 = /1 ® (-Z 2 ) <E)h® h, li = (tt, tt) T 4 = (-Zi) ® (-Z 2 ) <g> Z 3 ® 4, 
where, for example, the T 3 matrix is given explicitly by T 3 = diag[— 1, —1, —1, —1, 1, 1, 1, 1, —1, —1, —1, —1, 1, 1, 1, 1] 
and T4 = diag[— 1, —1, —1, —1, 1, 1, 1, 1, 1, 1, 1, 1, —1, —1, —1, —1]. The placement of the negative signs on the diagonal 
of Ti correspond to the 2 N ~ n blocks of size 2™ discussed in the Section and Appendix |[A"|). In this n = 2 example, 
there are 2 4 ~ 2 — 4 sections of size 2 2 = 4. 

The difference between T3 and the other tagging operations Ti is the location of the negative blocks on the diagonal 
of the T matrices that are associated with a matrix — Zi for each measurement angle Si = n on qubit i. Such an oracle 
is well behaved as defined by Bennet et. al [l6j]. Note that each tagging matrix Ti is, by construction, a tensor product 
of n 2 2 x 2 qubit matrices, and therefore a local, non-entangling operation. This differs from Grover's algorithm which 
can and in general will change the entanglement. 

We can now rewrite the fundamental quantum search equation Eq.© as, 

\y o ut,i) = H® n *T^ m ), (12) 

where T is the tagging operation chosen by the oracle for the desired output i e {1, . . . , 2™}. 

Interestingly, the tagging matrix Ti can be applied either to the right or to the left in Eq. (fl~2|l . In other words, we 
can apply the tag to the measurement operation (as typically done in MBQC) thus changing the iterators between 
iterations, or more interestingly, we can apply the tag to the initial input cluster state resource, prior to the mea- 
surement process. Applying the tag to the cluster state leaves the measurement sequence (pattern), or equivalently 
the MBQC search iterator, constant for any oracle chosen tag, regardless of the size N of the system, or the desired 
output of the search algorithm. The constant iterator is simply H® n acting on each column of qubits in turn. 

We now show the effect of acting with the tagging matrix Ti on the input state. Combining Eq.([T]) and Eq.© for 
the (n = 2) 4 element search we have, 

N 

= l^box) = n cz i!j \i>) = cz lt2 cz lt3 cz 2t4 cz 3A l[\+ i ). (13) 

i.j—nn i — 1 

Applying the Ti matrix we have, 

N 

Ti |* in ) = Ti \V box ) = Ti CZ h2 CZ h3 CZ 2A CZ 3A ]J \+ t ). (14) 
Ti and CZij are diagonal matrices, and hence commute we have, 

A' 

Ti |* m ) = Ti \V box ) = CZ h2 CZ h3 CZ 2A CZ 3A T l ]J (15) 

i=i 

Thus, we can say that the shape of the cluster state (it's pairwise connectivity), contained within the list of CZ gates, 
remains the same. The Ti operation only changes the state of the constituent qubits independently of the connectivity 
of the cluster state. We find that the effect of the Ti operation is to change several of signs of the initial cluster state 
from |+) states to |— ) states. There is a one to one correlation for a given measurement between those qubits, which 
in the standard implementation would be measured in the B(tt) basis, and the qubits which are rotated from |+) to 
I-)- 

As a further example, in the case of the (n = 3) 9 qubit MBQC search the measurement sequence (0, 0, tt ; 0, 7r, tt) 
on qubits (1, 2, 3 ; 4, 5, 6) corresponding to the fourth (of eight possible) output state (0, 0, 0, 1, 0, 0, 0, 0) T (see TableU). 
From Eq.([T]), the effective initial state for the search algorithm can be written as, 

N 

T 4 \il)i n u) = T 4 Y[ = |+i) <E> I+2) (8 |-s) <g> I+4) <8 |- 5 ) <8 he) ® I+7) <8 |+s) ® |+9>- (16) 
j=i 

Note the similarity of Eq. (TT51) to Eq. (l9bl) with — > Hi and |— —> Fj. These associations appears to hold true for 
any system size N and for any measurement patten (s\,s 2 ,s 3 ; S4,S5,sg). That is, any qubit measured in the B(n) 
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basis may equivalently be measured in the standard B(0) — |±) basis if the corresponding qubit in the initial cluster 
state is nipped from |+) — > |— ) with an overall phase which is dependent on the total number of Sj = 7r. 

For the (n — 3) N = n 2 = 9 qubit case, we explicitly display the form sequence of operations leading to MBQC 
search with a constant iterator. In the following we illustrate the specific example of tagging the fourth output state 
(i.e. again the U e i emen t4 measurement pattern) 

N 

Y^ouia) = U e i ement4 \^ in ) = (Hi <g> H 2 ® F 3 ) <g> (H 4 <g> F 5 <g> F 6 ) ® (H 7 ®H S ® H 9 ) JJ Cz iJ II ( 17a ) 

i,j—nn i — 1 

N 

= (Hi ® ff 2 ® H 3 ) ® (i? 4 ® ffs ® fle) ® (#7 ® # s ® #9) T 4 CZ,.,. (17b) 

i,j=nn i—1 

N 

= (H® n )® n n c/. j i , n ( i7c ) 

i,j—nn i—1 

n 

= H(ff® n ) fe H CZij |+i) ® |+ 2 > ® |- 3 ) <8 1+4) ® |- 5 > ® he) ® I+7) ® |+ 8 > ® l+ 9 ). (17d) 

fc=l i,j—nn 

Thus we have redefined the MBQC search in terms of a single tagging operation Tj which acts on a constant input 
state Ili=i l+i) followed by entangling operations CZij and a constant iterator defined as H® n acting on each of the 
n columns in turn (from left to right). In fully general terms the equations above become 

n N 

|«Wm) = l[(H® n )k II C Z *,3 T ml[\+l), (18) 
k—1 i,j—nn l — l 

where m is the desired output state in the decimal representation found by the tagging operation T m , N = n 2 is the 
total number of qubits n is the number of iterations as well as the size of one side of the square cluster-state and 
nn is the set of nearest neighbor in that cluster-state. In comparison, Grover's algorithm is quite different. Grover's 
oracle must repeatedly apply a unique tag to the state at each step so that the inversion about the mean operation 
can be applied correctly. Without this repeated tagging the algorithm would fail for converge to a unique answer. 
The MBQC algorithm can be performed with one, and only one, initial tagging operation to converge to a unique 
output state. 



VI. IS THERE AMPLITUDE AMPLIFICATION IN THE MBQC SEARCH ALGORITHM? 

We have now reformulated the MBQC search algorithm such that it utilizes a constant iterator, H® n applied to 
n qubit vertical columns, regardless of the required output state. Grover's iterator G = Ui nv Uf consists of two 
unitary operators: the (externally supplied) tagging operator U f chosen by the oracle, and the inversion about the 
mean operation, which were discussed above in section (IIV|) . The inversion about the mean operator Ui nv has entries 
Dij = 2/N G for i ^ j and (U inv ) u = -1 + 2/N G . This can be rewritten as U inv =2P' - I where P' = \ip in u) (Vw I 
is the projection matrix on to the unbiased state \ipinit) = ^2 X \x) with P[^ = 1/Nq for all The amplitude 
amplification takes advantage of the sign difference between targeted solution states after the tagging operation, and 
the re-adjustment of the amplitudes affected by the inversion about the mean operator. The Grover search algorithm 
requires repeated applications of the tagging operation Uf because Ui nv is Hermitian (Ui nv Ui nv = /), i.e. G 2 = I if 

In contrast, the MBQC search algorithm does not apply any amplitude amplification, nor does it repeatedly tag 
the state. As shown previously, the MBQC search oracle T only needs to tag the state once at the beginning of 
the algorithm. Recall that the MBQC iterator U described above implicitly contains projective measurements P 
similar to Eq. ©. The MBQC iterator is not an invertible operator for any size search algorithm. However, we can 
treat the MBQC search algorithm effectively implementing a unitary iterator on logical qubits, since the projective 
measurements performed on the columns of qubits in the n x n qubit grid leave the measured qubits in well defined 
states and the remaining logical qubits in a pure state. 

Brassard et. al generalized Grover's algorithm to a wider class of unitary amplitude amplification search problems 
[J, [13] ■ The general amplitude amplification iterator is defined by [l7l - [l9l | 



Q = AU A- 1 U f , 



(19) 
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Uf is the phase flip tagging operation discussed previously, Uq is the unitary operator that flips the sign of the 
computational basis state |0)® n while leaving all other computational basis states unchanged. The operator A is the 
generalized unitary oracle constrained only by the requirement that A|0)® n = \ipi n it)- Since the MBQC algorithm 
inherently requires projective measurements to proceed it can not implement amplitude amplification because it has no 
equivalent analogue to the inverse A -1 . We have confirmed this numerically by preforming incomplete MBQC searches, 
i.e. we preformed less iterations than the algorithm required for successful search completion, and examined the output. 
We found that the MBQC algorithm, which starts in the entangled cluster state \CS) = J\ i j =nn CZij\tpi n it) , does 
not alter the magnitude of any basis state amplitude, between the first and penultimate iteration (k = 1, . . . , n — 2). 
Thus no amplitude amplification takes place. 

Instead of amplitude amplification the MBQC search appears to be preforming a binary search operation. This 
empirical statement is backed by the computational time required to implement the MBQC algorithm. We find that 
the MBQC search algorithm requires exactly (n — I) iterations to find a solution with unit probability. Each MBQC 
iteration involves a simultaneous measurement of n qubits, for a total of n(n — 1) single-qubit measurements. If we 
include the final measurement to read out solution contained in the last column of qubits, the MBQC algorithm has 
applied n measurement iterations (measuring a total of n 2 qubits) to effect a search on 2 n elements. 

Note that any other smaller number of iterations (i.e. k < n — 1) will not yield a valid search solution, since the 
system remains in a state unitarily equivalent to the initial cluster state \CS). This result has been numerically 
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Number of Elements 2" 

FIG. 2: Iteration count vs number of searched elements, for classical and quantum unstructured search algorithms. The 
number of elements in an unsorted vector and the number of iterations required for three search algorithms. The classical 2™/2 
search (straight green line), the standard Grover's algorithm estimate of L(ti"/4)v / 2"J (red "stair-step" line) and the MBQC 
search algorithms log 2 2™ (black circles). (Colors Online). Note that the first point at which the MBQC algorithm becomes 
faster is 59 elements, which requires a 6x6 qubit square grid. The last point at which Grover's is as faster as the MBQC search 
is an 80 element search on a 7x7 qubit square grid. 

tested for searches of 2™ elements with n = 2, . . . , 9. The (n = 9) 2 9 = 512 element MBQC search requires 81 total 
qubits, 72 single qubit measurements and up to 72 post rotations. Such a large cluster state is currently physically 
un-realizable, but our numerical results show that such a cluster state takes only n — 1 = 8 applications of the MBQC 
iterator compared to (7r/4)v2™ s» 17 applications of Grover's iterator. This result is reflected in Fig. ^ where the 
abscissa has been truncated to maintain clarity for the smaller element sizes which are currently achievable. 

This suggests that the MBQC search is not amplitude amplification, rather a l-to-2™ mapping, and that the oracle 
is simply selecting a given path. As each single-qubit measurement in the MBQC search has only 2 outcomes, the 
MBQC search algorithm can be thought of as a binary tree structure. The oracle directs the search through the binary 
tree (from root to leaf) based on the oracles application of the tagging operation. Each individual measurement in 
either the or ir rotated basis (i.e. B(0) or B(?r) from Eq.(|8|)) selects a branch of the binary tree to move down. As 
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stated above, some measurements are not necessary if one ignores or corrects for overall phase effects. In fact all 2™ 
outcomes for the MBQC search algorithm can be found with a minimum of exactly n variable measurements. This is 
suggestive of the binary structure, since the total number of possible unique outcomes is 2™. 



VII. ENTANGLEMENT DYNAMICS 
A. Relative entropy of entanglement 

As another quantitative means of comparison between the MBQC algorithm and Grover's algorithm, we calculated 
the relative entropy of entanglement for the logical three qubit and four qubit searches. For Grover's algorithm this 
is relatively straightforward. We need only record the system state vector \ip) at any given stage of the calculation, 
form the density matrix p — \^p){ip\, and calculate the relative entropy relative to the closest separable density matrix 
a . The MBQC calculation is slightly more involved as one must identify the current n qubit logical pure state within 
the n x n qubit cluster state during the computations, after each projective measurement. 

We calculate the three qubit and four qubit relative entropy of entanglement for each step of the respective searches 
numerically via a numerical optimization routine [20| . The relative entropy of entanglement is based on the Von 
Neumann entropy, 

S(p) = -Tr(plnp). (20) 

We use the relative entropy of entanglement because we desire a measure of the distinguishably of our two search 
states compared to the set of all separable states. The relative entropy of entanglement is defined as [2Ct |2l|. 

E q (p) = m% 6 s, - Tr(plnp - phxx), (21a) 

= mhx^sM -lna|V) (21b) 

where a is an element of the set of all g-separable states S q , namely states that are separable for q partitions of 
the n-qubit system, and p is the density matrix representing the logical state of the system [22| . In Eq. (|21b[) we 
have used the fact that p = is a pure state density matrix. For E n , a is n-qubit separable. E n quantifies the 

amount of entanglement of any kind contained in the global system, i.e. it is non-vanishing even for states showing 
ent ang lement just between two subsystems (2ll . [23| . For E2, which quantifies genuine multipartite entanglement 
[HI, HH, the minimization of a takes place over all possible bipartite partitions. 

For the n = 3 qubit case, there are only three possibilities, (i) E\ = (since a is a pure 3-qubit state), (ii) E2 in 
which a is bipartite separable and (ii) E3 in which a is 3-qubit separable: We define the set of separable states a 



numerically as, 

^=2) = X>(l&,i,a> I'M) «&,1,2|<<M)> ( 22a ) 



i=l 



N. 



(9=3) 



a 



$> (l<MI<MI<M) (W> M |(<M(<M), (22b) 



i=i 



where pi is a probability with J^Pi = 1j an d No = 2™ is the number of distinct search elements. In Eq. (|22a[) 1^,1,2) 
represents a general two-qubit state on qubits 1 and 2 and |^s) represents a general single-qubit state on qubit 3. 
The other possible grouping of qubit labels, i.e. {(1,3), 2} and {(2, 3),1} differ from p^zS by local unitary swap 
operations of the form U± <S> U2 <8> U3, and hence are effectively searched over in the minimization procedure due to 
the symmetry of the quantum search problem. In Eq. (|22b|) \<f>j t i) is a general single-qubit state on qubit j e {1, 2, 3}. 
Note that for n qubits the set of all density matrices is dimension Nq x Nq, and therefore any density matrix can be 
written as a linear combination of Nq basis matrices (e.g. (Mij) a ^ = 5i <a 8j t b, for a,b G {1, . . . , Nq}). 

For the n = 4 case, E3 is a straightforward generalization of Eq. (|22a[) with an extra single-qubit partition on qubit 
4, and E4 is a straightforward generalization of Eq. (|22b[) now involving four single-qubits. The relative entropy E2 
now involves two types of bipartite partitions of the form 

a (q=2) = ^Pi(\4>ia^)\^A))(( < t>hh!i\(<l>i,3A\), (23a) 
i=l 

a (S~3) } = E^(l^^)l^,i»(^l,2,3,i|^4,i). (23b) 



i=l 
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Both forms of separable density matrices in Eq. (I23a[) and Eq. ()23b[) must be searched over in the minimization 
procedure. We calculate the relative entropy of entanglement value by numerically minimizing over the set of separable 
states a and take advantage of the fact that the relative entropy is convex thus a local search results in a global 
minimum. We have preformed this calculation for each step of the two search algorithms for both the n=3 qubit 8 
element search and the n=4 qubit 16 element search. 



B. The pure state system density matrix 



The use of a pure state system density matrix p in the standard Grover's algorithm this is straightforward as we are 
dealing with a pure state vector of size 2™ at each iteration. As stated above for the MBQC search the use of a pure 
state p is slightly more complicated and requires some interpretation regarding the logical n qubit state in the space 
of the larger cluster state. From Raussendorf, Browne and Briegel [HJ we learn that the MBQC algorithm can be 
broken into a series of projective measurements and entangling operations. This can be seen in the following relation, 



{P 2 D 2 ){P 1 D 1 ) = [P 2 P 1 ){D 2 D 1 ) 



(24) 



Here {Pi} is a set of n qubit projective measurements, thus removing a column of qubits, and {Di} is a set of 
entangling operations 24| acting on the logical state and adding new qubits to the cluster state. The key point to 
note is that the ZVs commute with the Pfc's. This fact allows for the piecewise construction and computation on 
cluster states. In other words, one can implement the MBQC search algorithm (ignoring whether or not it is the same 
as Grover's algorithm) while never requiring the presence more than 2n qubits at the same time; the additional qubits 
can be entangled "on the back end" and operated upon at a later time. Stated slightly differently, the feed forward 
operations in the MBQC implementation never effect more than nearest neighbor qubits. We used this feature in our 
calculation of larger cluster state such as the 81 qubit cluster state discussed in section (TVip . 
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FIG. 3: Schematic of piecewise cluster state implementation. Cluster state map and effective pure state vector in 
MBQC algorithm (see text for description). 



In the (n = 3) 2 3 element MBQC search implemented on a 3 x 3 grid of qubits illustrated in Fig. ([3]), the hrst 
three qubits 1,2,3 are initialized into the |+) state yielding the unbiased state \tpmit) (the same state utilized in the 
3-qubit Grover search algorithm at iteration k — 0). Subsequently, the qubits are CZ-entangled with each other, 
yielding the pure state |-0o) = C Zi y2 C Z 2i3 \tfj in it) = \C3). \i/}q) is the initial system pure state of the MBQC search 
algorithm at iteration k = 0, and is the three qubit linear cluster state |C3). We expect this state to be highly 
entangled, which can be seen in Fig. (Hk), where we plot the relative entropy E 2 for both Grover's algorithm (red) and 
the MBQC algorithm (black) for 3 qubits, versus the iteration number k. The MBQC search iterator then begins its 
first iteration (fc = 1) by entangling |V>o) with qubits 4,5,6, which are also begin in the |+) state, and are themselves 
CZ-entangled with each other, i.e. D\. This creates a 2 x n cluster state, similar to the butterfly state [25]. The 
qubits 1,2,3 are then projectively measured, i.e. P±. This leaves a three qubit state on 4,5,6 which we denote as the 
pure state \ipi). (Recall that without loss of generality we can consider qubits i = (1, 2, 3) to be projected into their 
respective |+ s ,i) state, with any required Pauli corrections taking place in post-processing). This is now the three 
qubit system pure state after one iteration (k = 1) of the algorithm. In the next iteration (k — 2) the state \ip\) 
is entangled with three new qubits 7,8,9, identically to the previous step. For clarity we call this unitary operation 
D 2 . Qubits 4,5,6 are then measured, P 2 , and the logical pure state ^2) is contained in qubits 7,8,9. We note that 
the projective operations Pi are simply the three qubit projective iterators discussed in section (jVj, and the Di are 
the piecewise construction, i.e. entanglement, of the cluster state illustrated in Fig. ([3]). The qubits 7,8,9 need not 
exist when qubits 1,2,3 are measured since the quantum information is entirely stored in 4,5,6. The consequence is 
that no more than 2n qubits are ever needed at any given time in the algorithm. This has particular relevance to 
cluster states constructed from "flying" qubits. The significance is that no more than 2n photons need ever be present 
simultaneously. This constitutes a significant advantage for the photon-based MBQC architectures , especially since 
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the creation of large photonic cluster states is currently limited by the difficulty inherent in producing larger number 
of on-demand single photons. 

The above piecewise constructive procedure also yields a clear sequence of three and four qubit logical pure states, 
HV'o); l^i); 1^2), •■•} that we can compare directly with corresponding pure states produced in the standard Grover's 
algorithm. We have carried out this comparison using the relative entropy of entanglement Ei as defined above. The 
results are shown in Fig. (j4j. Here the number of iterations has been increased beyond the ideal number of steps for 
successful search completion for the 8 element and the 16 element searches in order to reveal more behavior of the 
algorithms. 
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FIG. 4: Relative Entropy of entanglement E2 vs. iteration number k for (a) the three and (b) the four qubit logical 
states of both quantum search algorithms. Red squares, the relative entropy of entanglement of Grover's Algorithm. Black 
circles, the relative entropy of entanglement of the MBQC search algorithm. Note that while the red curve in figure (a) appears 
symmetric, it is not, nor are the states with similar E2 identical. 



C. Relative entropy analysis for the Grover search algorithm 



We expect the entanglement of the standard Grover's's algorithm to be E% — for the input state, an n-qubit 
separable product state, and then increase to a maximum after L(ti"/4)v / 2™J steps (where [x\ is the nearest integer 
less or equal to x) before returning close to zero again on the solution state (approximately an n-qubit separable 
product state). This was shown by M. Rossi et. al [21(. This is confirmed by our results for the first 2 steps (Grover 
iterations of U[G) in Fig. and the first 3 steps of Fig. <m>) (red). Note the graph for the the standard Grover's 
algorithm doesn't return to E2 — at the ideal terminal iteration (k = 2 for the 8 element search, and k = 3 for the 
16 element search) where the solution is achieved, because the standard algorithm is not perfect at this step, and we 
are not performing any projective measurements on the output states. Iterating past the solution step, the spike in 
the entanglement at step k = 3 in Fig. (j4^,) and step k = 4 in Fig. (j4j)) is due to the imperfect nature of Grover's 
algorithm. It is interesting that the most entangled state is never utilized in Grover's original algorithm. G. L. Long 
proposed a modified Grover's algorithm that is almost perfect (i.e terminates on a solution that is almost exactly an 
n qubit separable state) for larger size systems @, that should eliminate this spike. 

The general behavior of Grover's algorithm can be understood from it's periodic behavior with produces the state 
IV'fc) = Uq \ipinit) after the fc-th application of the Grover iterate XJq- As is well known [l7l - [T9j . the unbiased state can 
decomposed into the two-state basis \ipinit) = sm9\ip goo d) +cos6\ipbad}- Here, \ip goo d) — \xq) is the normalized solution 
state (chosen to be 1 0) 1 |1) 2 |0) 3 for the graph in Fig. (@ji)), and \ipbad) = v/jv ^_ 1 J2 X jtx witn sin ^ = is tlic 
superposition of all the non-solution states. After k iterations of Ug, the resulting state is given by \ipk) = U G \ipi n it) 



13 



= sin(2fc + l)0\ipgood) + cos(2fc + l)0\ipbad), with the optimal number of steps for a solution given by the condition 
that (2k + 1)8 = n/2 yielding the above estimate kcrover ~ L( 7r /4)v / 2"J • In general, a Grover search solution will 
be obtained under the condition (2k + 1)8 = (2m + l)ir/2, for any integer m. In Fig. (QJt) the second occurrence of 
E 2 ~ in the graph for Grover's algorithm occurs for k = 6 and k = 8 for the 3-qubit and 4-qubit search respectively, 
yielding estimates of the integers m (utilizing 8 rj {l/VS, 1/v16}) of m w {0.987,0.902}, both roughly close to 
m = 1. The peak (spike) £2 ~ 1 occurs at fe = 3 and k = 4 for the 3-qubit and 4-qubit search respectively, with 
(2k + 1)8 = (0.7167T, 0.7887t). These values are both roughly 37r/4 associated with the maximally entangled state 
\ipE 2 =i) = (\^l J good) — l^bad))/^, which has the same magnitude for all amplitudes (as in the unbiased state l^mit)), 
but with an important relative minus sign between the solution and non-solution states. Finally, the secondary 
maxima of E2 for Grover's algorithm for the 3-qubit case occur for iterations k = 1 and k = 5. At these steps, all 
computational basis state amplitudes have the same sign, and all have equal magnitude except for the solution state 
component whose magnitude is significantly larger (but whose square is not close enough to unity for the whole state 
to be approximately \xq)). Similar behavior occurs for the secondary E2 maxima in the 4-qubit case. 



D. Relative entropy analysis for the MBQC search algorithm 

We expect the MBQC state, which starts in the 3-qubit case in a highly entangled state l^o) = |C3) = 
CZi j 2Ci?2.3|+)i|+)2|+)3, to decrease in entanglement as the projective measurements are made. Here we'll write 
the state vector at the A:-th iteration \ip k ) in the |±) representation so that |^q) = [1 ; 0; 1> 0) 0, 1, 0, — 1]/V8 with the 
first component corresponding to |+) 2 |+)3 and the last component corresponding to |— )i |— ) 2 |— }a. In relation 
to Grover's algorithm which begins in unbiased state \ipo) — \ipinit), the MBQC puts all the entanglement in the initial 
state at iteration k = with E 2 = 1 as shown the black curves in Fig. (Q| . For the 3-qubit case shown in Fig. (with 
tag T3 selected, i.e. a constant iterator of (7r,7r,7r), see the last row of Table llVl in Appendix |C|) . the relative entropy 
of the state \ip\) = 1, 1, 1, 1, —1, 1, 1, 1]/V8 is unchanged after the first MBQC iteration and remains at E2 = 1. 
At the second (n — 1 = 2) iteration the state \tp2) is transformed into the solution state \xo) = [0, 0, 1, 0, 0, 0, 0, Oj/v 7 ^ 
= l+i) |— 2) I+3) = H® 3 |0i) 1 12) IO3) with E 2 = 0. Continuing to iterate past the successfully completed search, the 
state after the third iteration \ip 3 ) = i[— 1, 1, — 1, 1, 1, — 1, 1, 1, — l]/s/8 = F® 3 ]^) again has E 2 = 0. At the fourth 
iteration, we calculate ^4) = [—1, 0, 1, 0, 0, —1, 0, — l]/v^ with E 2 = 1, very similar in structure to the initial cluster 
state |C3). Additionally, the two iterations k = 4, 5 produce states similar in structure and relative entropy E2 = 1 
as in iterations k = 0,1. In fact, the states at k = {0,1,4,5} are local unitary equivalent to each other, i.e. differ 
on by local unitary transformations of the form U\ ® U2 <8> U3. Finally, iteration k — 6 returns again to a relative 
entropy E 2 = with output state \tp 6 ) = [0,0,0,0,0,0,0,1]/^ = |-)i|-) 2 h>3 = H® 3 |l)i |1) 2 |1) 3 - The state 
is the computational basis solution state that would have arisen at iteration k = 2 if one had initially implemented 
the tagging operation Tg (see the third row of Table HVl in Appendix [C]). The states at k = {2, 3, 6} are local unitary 
equivalent to each other. Qualitatively similar behavior is exhibited in the 4-qubit case as show in Fig. (J3t>) - 

The above cyclic behavior of the "solution sets" of the MQBC algorithm is evident from Table HVl for the 3-qubit 
case and Table IVl for the 4-qubit case in Appendix ICl In these tables we repeat the entangling/rotation/projection 
operations many steps beyond iteration k = n — 1 where the solution associated with the chosen tag first appears, i.e. 
the state \x = i), in the decimal representation, associated with tag for i g {1, . . . , 2™ — 1}. For a given tag Tj, the 
solution either (i) stays the same (modulo overall phase) or (ii) periodically cycles between two states (modulo factors 
of ±1 and/or ±z) every n + 1 iterations. Note that the "orbits" for different tags may cross on these solutions. That 
is, a given solution state \x) (e.g. x — (2, 4, 5, 7) in Table (|IV|I ) may appear in the "orbits" of more than one tag due 
to the differences in the measurement bases involved in the MBQC iterators for each chosen tag. 

The MBQC states in Fig. (|4^l) and Fig. (|4Jd) with E 2 — 1 can be understood as follows. For the 3-qubit case at k = 
the initial state is the cluster state |C3) = (|0 X ) |C2 2)3 ) + |li) \C2 2 J)/V2 with |C2 2 , 3 ) = (I+2) |0 3 ) + |- 2 ) |1 3 >)/V^ 
and |C2' 2 3 ) = (|— 2) | O3) + I+2) 1 13) )/ v^- A representative bi-separable state yielding E2 = 1 is given by ^[^2) = 
l/2(|0i, C22,3)(0i, C22,3| + 1 1 x , C2' 2 3 )(li, C2' 2 3 |), which is the sum of the projectors associated with each of the two 
terms in |C3). This state is bi-separable in the qubit partition {1,(2,3)}, while exhibiting maximal entanglement 
on the (2,3) sub-partitions embodied in IC22.3) and |C2' 23 ). All other bi-separable states with E 2 = 1 are local 

unitary equivalent to cr^Z 2 y The first measurement does not change the value of E 2 so that the state at k = 1 
is again local unitary equivalent to a9^~S. Similar behavior occurs in the 4-qubit case. At k = the cluster 
state is given by |C4) = (|0i) |+ 2 ) |C2 3)4 ) + |li) |- 2 ) \C2% 4 ))/y/2 with |C2 3 , 4 ) with the same form as before and 
|C2 3 ' 4 ) = (I+3) I O4) — |— 3) 1 14) )/ a/2- A representative bi-separable state yielding E 2 = 1 is given by o - /™^) = 
l/2(|0i,+2,C2 3i 4)(0i,+2,C2 3 ,4| + \l 1 ,- 2 ,C2'l 4 )(li,-2,C2'l 4: \). Again, all other bi-separable states with E 2 = 1 are 
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local unitary equivalent to the states at k — 1 and k = 2. 



VIII. REMAINING QUESTIONS 

We have examined a number of aspects of the MBQC search algorithm, however, several questions remain unan- 
swered. How can a search be preformed that has more than one tagged output in the MBQC algorithm? How does 
this affect the speed of the algorithm? Can the idea of critical damping be applied to such multi-element MBQC 
algorithm? Can the database search algorithm on "itemized lists" encoded in tensor products states as described at 
the end of Section IIV Al 15] be implemented in a MBQC algorithm? Can the MBQC algorithm be implemented on 
a blind cluster state [26|, 127J with similar speed? These questions will be explored in future work. 

IX. CONCLUSION 

Grover's Algorithm and the MBQC search algorithm (that has been called Grover's algorithm in the literature [Toj ) 
are qualitatively similar in form but significantly different in structure and operation. Both search algorithms appear to 
be oracular search methods, mapping a chosen tagging operation, performed by the oracle, with a corresponding unique 
search solution. However, the behavior of the oracle in each search algorithm is quite different; in Grover's algorithm 
the tagging is an entangling operation and is performed at each iteration, while in the MBQC search algorithm the 
tagging operation is non-entangling and can be performed on the unbiased initial state, before the application of 
the entangling CZ gates used to create the initial cluster state. The iterators of the two search algorithm are also 
structurally very different. In Grover's search algorithm, the iterator preforms amplitude amplification, while in the 
MBQC search algorithm the iterator appears to be a selection a of single "route" or "path" in a l-to-2" mapping of an 
initial state the set of all possible solution states. This is most readily exhibited through the entanglement dynamics 
presented in this work of the n-qubit logical state at each iteration. We have shown explicitly that the iterators in 
the MBQC search algorithm do not perform amplitude amplification; the projective measurements toggles between 
states of the same entanglement class, either similar to the entangled initial state, or similar to the n-qubit separable 
solution state. Significantly, the speed of the two search algorithm are very different, with the MBQC search being 
the faster of the two. That is, for a 2 n element search, Grover's search algorithm requires 0{2 n ' 2 ) applications of the 
Grover iteration acting on n-qubits. For the MBQC algorithm acting on n 2 qubits, only n — 1, i.e. 0(n) applications 
of the MBQC iterator are required, acting on ro-qubits at each iteration. From the above analysis and observations we 
draw the conclusion that MBQC search algorithm is not Grover's algorithm, and should be studied in its own right. 
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Appendix A: Pauli corrections and the projective measurement P 



In Section (JTTTJ) we state that the output matrix \ip ut) — U\ipo) is a vector of size 2 N (see Eq. ([5])). We then 
claim that without loss of generality we can consider the logical output state to always be contained within the 
first 2™ elements of the output vector, and ignore the remaining 2 N ~ n components. This statement follows from 
the non-deterministic projective nature of quantum mechanical measurements. Whenever we attempt to make a 
measurement on a superposition of two states (say the |+) and |— ) states) we can obtain either output. In the MBQC 
implementation one makes a partially ordered series of single-qubit measurements on the |+ s ) state (|+o) = |+) and 
= \— ))• For each single-qubit measurement there is equal probability that a particular outcome will yield the 
"incorrect" |— s ) measurement, (|— o) = | — ) and |— „) = |+)). This does not present a problem since the incorrect 
measurement merely adds a trivial phase and Pauli flip X to the desired target state. Eq. (IA1|) shows the effect 
described in the main body for the case of the n=2, 4-element, search, where measurements will be performed on 
qubits 1 and 2 as in Fig. [U with the output state projected on qubit 3 and 4. Equation (IA1|1 show the entire 4-qubit 
state prior to the projection operation. Here we use the unitary matrix (tagging and iteration) that yields the 4th 
element (0,0,0, 1) T on qubits 3 and 4 (analogous to Eq.((9b|). 






1 



\1pout) = V elementiWo) 




_0_ 

-1 




V o y 




(Al) 



Notice that the desired logical output is indeed contained with the first 2™ elements of \tp ut) (the "first section"). The 
remaining 2 n terms contain n — 1 copies (sections) of the targeted state \tptarget) (on qubits 3 and 4), each of which 
has undergone a unique set of Pauli bit flips and overall phase manipulations corresponding to those single-qubit 
measurements (on qubit 1 and 2) that were projected into the incorrect states |— s ) . In effect we are calculating every 
possible outcome of the n(n — 1) measurements simultaneously, and therefore the state before the projection is an 
equal superposition of all these possible outcomes. An experiment would produce only one of these four output states 
with equal probability. Since we know the outcome of the measurements (i.e. we observe the measurement record), 
we always know which particular qubits were projected into the state |— s ). Thus, we can always construct the Pauli 
bit flip and overall phase operations shown in Eq. (|A1|) . to correct (in post processing) the output such that we end 
up with exactly \iptarget) 011 qubits 3 and 4. For operational simplicity and computational speed, we truncate our 
matrices U to the first 2" rows of the full 2 N x 2 N unitary matrix, thus containing one copy of the "correct" (i.e. not 
Pauli flipped or phase shifted) logical output. This is implemented numerically by simple truncation of the matrix U 
and justified logically by stating that we can assume, without loss of generality (since we can perform the necessary 
Pauli corrections in post processing), the measured qubits are always projected into the correct |+ s ) states. This is 
equivalent to applying the operator P = |+«)(+ a | (identical to Eq.©) on each each measured qubit. 



Appendix B: Measurement setting for the MBQC search 



As illustrative examples we provide measurement settings for the MBQC search for the (n = 3) 8 and (n = 4) 16 
element searches. Table |TIJ shows the 8 element search similar to the main body, but additionally also shows the 
setting which provide constant iterators. Observe that in both sets of measurements there are only 3 free variable 
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measurements. Note the the (trivial) overall phase factors which are ignored in this table are not the same between 
the two measurements settings. This demonstrates part of the n-fold degeneracy that was mentioned in the section 
on the tagging operation. 



TABLE II: Measurement sequences & resulting output ignoring phase factors for 2 3 elements 



Constant Iterators (T\ ; Zjj) 


Variable Iterators 


Element i 


Output State \&out,i) 


(0,0,0 


0,0,0) 


(0,0,0 


0,0,0) 


1 


(1,0,0,0,0,0,0,0) 


(0,7T,7T 


0,7T,7r) 


(0,0,7T 


0,7T,0) 


2 


(0,1,0,0,0,0,0,0) 


(7r,7r,7r 


7r,7r,7r) 


(0,0,0 


0,0,7T) 


3 


(0,0,1,0,0,0,0,0) 


(tt,0,0 


7T,0,0) 


(0,0,7T 


0,7T,7r) 


4 


(0,0,0,1,0,0,0,0) 


(7T,7T,0 


7T,7T,0) 


(0,0,7T 


0,0,0) 


5 


(0,0,0,0,1,0,0,0) 


(7r,0,7T 


7T,0,7r) 


(0,0,0 


0,7T,0) 


6 


(0,0,0,0,0,1,0,0) 


(0,0,77 


0,0,7T) 


(0,0,7T 


0,0,7T) 


7 


(0,0,0,0,0,0,1,0) 


(0,7T,0 


0,7T,0) 


(0,0,0 


0,7T,7r) 


8 


(0,0,0,0,0,0,0,1) 



Table fjlllf) shows the 16 element search and the constant iterator measurement set. Here we are able to show the 
phase factors for these measurements. 

TABLE III: Measurement sequences & resulting output including phase factors for 2 4 elements 



Constant Iterators (T\ \Ti ; Ta) 


Element i 


Output State \^ ut,i) 


(0,0,0,0 


; 0,0,0,0 ; 


0,0,0,0) 


1 


(1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0) 


(7r,7r,7r,0 


7T,7r,7r,0 


7r,7r,7r,0) 


2 


(0,-i,0,0,0,0,0,0,0,0,0,0,0,0,0,0) 


(i",0,vr,7r 


7r,0,7T,7T 


7r,0,7r,7r) 


3 


(0,0,-i,0,0, 0,0,0,0,0,0,0,0,0,0,0) 


(0,7T,0,7T 


0,7T,0,7r 


0,7T,0,7r) 


4 


(0,0,0,-1,0,0,0,0,0,0,0,0,0,0,0,0) 


(7T,7T,0,-7r 


7r,7r,0,7r 


7r,7r,0,7r) 


5 


(0,0,0,0,-i,0,0,0,0,0,0,0,0,0,0,0) 


(0,0,7T,7T 


0,0,7T,7r 


0,0,7T,7T) 


6 


(0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0) 


(0,7T,7r,0 


0,7T,7r,0 


0,7T,7r,0) 


7 


(0,0,0,0,0,0,-1,0,0,0,0,0,0,0,0,0) 


(tt,0,0,0 


7T,0,0,0 


7T,0,0,0) 


8 


(0,0,0,0,0,0,0,-i,0,0,0,0,0,0,0,0) 


(0,7r,7T,7r 


0,7r,7r,7r 


0,7T,7r,7r) 


9 


(0,0,0,0, 0,0,0,0,-i,0,0,0,0,0,0,0) 


(7T,0,0,7T 


7T,0,0,7T 


7T,0,0,7r) 


10 


(0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0) 


(7T,7T,0,0 


7T,7r,0,0 


7T,7r,0,0) 


11 


(0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0) 


(0,0,7T,0 


0,0,7T,0 


0,0,7T,0) 


12 


(0,0,0,0, 0,0,0,0, 0,0,0,-i,0,0,0,0) 


(7r,0,7T,0 


7r,0,7T,0 


7r,0,7T,0) 


13 


(0,0,0,0,0,0,0,0,0,0,0,0,-1,0,0,0) 


(0,7T,0,0 


0,7T,0,0 


0,7T,0,0) 


14 


(0,0,0,0, 0,0,0,0, 0,0,0,0,0,4,0,0) 


(0,0,0,7T 


0,0,0,7T 


0,0,0,tt) 


15 


(0,0,0,0, 0,0,0,0, 0,0,0,0,0,0,-i,0) 


(7r,7r,7r,7r ; 




; 7r,7r,7r,7r) 


16 


(0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-1) 



Appendix C: Cyclical behavior of the constant iterator MBQC search algorithm 

We find that the MBQC search is also cyclical in nature when a constant iterator is applied. For a given tag, 
Grover's algorithm cycles between the desired target state, an n-qubit separable state, and a different solution state 
which approximates an rt-qubit separable state. In contrast, the behavior of the MBQC search is more complicated. 
For a given tag, the desired target state will be produced at iteration k = n — 1. Afterwards, for some tags, the desired 
target state will be produced every n + 1 iterations. However for other tags, the states cycle every n + 1 iterations 
between the target state and a second distinct valid output state that is exactly an n-qubit separable state. Note that 
the desired target state repeats with period 2(n+ 1) regardless of the chosen tag. Table (|IV[) shows the chosen iterator 
and the resulting tagged element from 1 to 2 3 = 8 for the case n = 3, including overall phases at the output of each 
cycle. Table (jVj) shows the corresponding results for the case n — 4. The output states in Table (|IVp are similar to 
Fig. (J4^i) at iterations k = 2 and k — 6, and beyond (similarly for Table ([Vj) and the output states at iterations k = 3 
and k = 8 in Fig.(j4j})). Thus, each column requires a different size cluster state or a different number of piece wise 
operations Pi and Di similar to Eq. (|2~4"l) . For the n — 3 case, the iterators (0,0,0) and (tt, 0,tt), corresponding to the 
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TABLE IV: The tagged element and the number of iterations at which it is reached for the 2 3 element search. 



Constant Iterator X 


Iteration 2 


Iteration 6 


Iteration 10 


Iteration 14 


Iteration 18 


(0,0,0) 


1 


1 


1 


1 


1 


(0,0,tt) 


-7 


2 


7 


-2 


-7 


(0,7T,0) 


8 


-3 


8 


-3 


8 


(0,7T,7r) 


-2 


-4 


2 


4 


-2 


(tt,0,0) 


-4 


5 


4 


-5 


-4 


(7T,0,7r) 


6 


6 


6 


6 


6 


(7T,7r,0) 


-5 


-7 


5 


5 


-5 


(7T,7r,7r) 


3 


-8 


3 


-8 


3 



target states 1 and 6 respectively, give output elements that are always identical. The remaining tags cycle between 
between two n-qubit separable target states. The precise cause of this cycling behavior is a topic warranting further 
investigation. 

TABLE V: The tagged element and the number of iterations at which it is reached for the 2 4 element search. 



Constant Iterators X 


Iteration 3 


Iteration 8 


Iteration 13 


Iteration 18 


Iteration 23 


Iteration 28 


(0,0,0,0) 


1 


1 


1 


1 


1 


1 


(7r,7r,7r,0) 


-15i 


2 


15i 


-2 


-15i 


2 


(7r,0,7r,7r) 


-12i 


3 


12i 


-3 


-12i 


3 


(0,7T,0,7r) 


6 


4 


6 


4 


6 


4 


(7r,7r,0,7r) 


-14i 


5 


14i 


-5 


-14i 


5 


(0,0,7T,7r) 


-4 


6 


-4 


6 


-4 


6 


(0,7T,7T,0) 


-7 


7 


-7 


7 


-7 


7 


(tt,0,0,0) 


-9i 


8 


9i 


-8 


-9i 


8 


(0,7r,7r,7r) 


-8i 


9 


8i 


-9 


-8i 


9 


(7T,0,0,7r) 


10 


10 


10 


10 


10 


10 


(7T,7r,0,0) 


-13 


11 


-13 


-11 


-13 


11 


(0,0,7T,0) 


3i 


12 


-3i 


-12 


3i 


12 


(7T,0,7T,0) 


11 


13 


11 


13 


11 


13 


(0,7T,0,0) 


5i 


14 


-5i 


-14 


5i 


14 


(0,0,0,tt) 


-2i 


15 


2i 


-15 


-2i 


15 


(7r,7r,7r,7r) 


-16 


16 


-16 


16 


-16 


16 



